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The conventional Hamming distance measurement captures only the short-time dynamics of the 
displacement between the uncorrelated random configurations. The minimum difference technique 
introduced by Tirnakli and Lyra [Int. J. Mod. Phys. C 14, 805 (2003)] is used to study the 
QO ' short-time and long-time dynamics of the two distinct random configurations of the isotropic and 

anisotropic Bak-Sneppen models on a square lattice. Similar to 1-dimensional case, the time evo- 
lution of the displacement is intermittent. The scaling behavior of the jump activity rate and 
CNj ' waiting time distribution reveal the absence of typical spatial-temporal scales in the mechanism of 

» I | displacement jumps used to quantify the convergence dynamics. 
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^ 1 I. INTRODUCTION 

o 



The term self-organized criticality (SOC) was first introduced by Bak et al. in 1987 In their well-known paper, 
it was argued that the dynamics which gives rise to the power-law correlations seen in the non-equilibrium steady 
states must not involve any fine-tuning of parameters. Namely, the systems under their natural evolution are driven to 
a statistically stationary state where long range spatio-temporal fluctuations are seen similar to those in equilibrium 
critical phenomena. 

The Bak-Sneppen (BS) model, which models the biological evolution of an ecology of interacting species Q is one 
of the simplest model that exhibits SOC. Besides its simplicity, this prototype of SOC shows a quite rich variety of 
properties from the critical phenomena. In the BS model, an ecosystem consisting of L d species on a d-dimensional 
■ lattice is characterized by random fitness values, each of which represents the adaptation of each species to the 
£h ' environment. 

In this paper, we consider only 2-dimensional case. The dynamics of BS model on a 2-dimensional lattice of an 
edge size L has simple rules. The initial state of the system on a square lattice is characterized by N = L x L fitness 
values (random numbers) fi,j, where i = 1 • • • L and j = 1 • • • L, uniformly distributed between and 1. These fitness 
values are assigned to each site (i, j) of the lattice with periodic boundary conditions. The usual dynamics of the 
system is eventually achieved by localizing the lattice site with the minimum fitness / m j n and assigning new random 
£^ i numbers to that site and its first nearest neighbors. 

The model described above can be called as an isotropic BS model since the interaction between the current 
i minimum and its first nearest neighbors is the same in both directions. In other words, if one observes the minimum 
fitness value of the system at time step t as f min = faj then the fitness values fi-i,j, fi,j-i, fij, fi+i,j an d fi,j+i 
will be updated at that time step. This means that the possibility for the minimum to jump to one of its left or right 
nearest neighbors and to one of its up or down nearest neighbors in the next time step is simply the same, with a 
resulting isotropic avalanche of events. On the other hand, one can easily consider alternative updating rules. One 
• • i possible alternative could be to update at each time step fi- Ut j, fi.j-i, f%,j, fi+b,j and fi,j+ r (where u, £, b and r 
are arbitrary positive integers taken from the interval [1, L]). Let us consider, for example, the case of fi-2,j, fi,j-i, 
fi,j, fi+i,j an d fi,j+2- The system now has an inherent bias to the up and right. This means that we would expect 
a preferred direction for an avalanche to propagate. Therefore, this model is called as the anisotropic BS model. In 
principle, several types of anisotropy can be introduced by changing the values of u, I, b and r, provided that £ r 
and u ^ b (since £ = r and u = b represents the isotropic BS model). The maximal anisotropic cases are defined by 
I = (£ = 1), r = 1 (r = 0) and b = (b = 1), u = 1 (u — 0) whereas other definitions arc considered as intermediate 
anisotropics. In this work, we use the maximal anisotropic case in all simulations since the convergence is faster than 
the intermediate anisotropy choices 0]. 

After some transient time which depends on the size of the system the isotropic and anisotropic models achieve 
a statistically stationary state (i.e., self-organized critical state) in which the density of fitness values is uniformly 
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distributed on [f c , 1] and vanishes on [0, f c ], where the critical threshold value f c depends on the lattice dimension. 
On the 2-dimensional lattice considered here, f c ~ 0.328 3 for the isotropic BS model and f c ~ 0.439 [5] for the 
anisotropic BS model. The difference between the critical threshold values of isotropic and anisotropic BS models 
comes from the change in the rates of the spreading out of avalanches. Once the stationary state is achieved, the 
temporal and spatial correlation functions are power-law in both models signifying the existence of a critical state with 
no characteristic length or time scales (scale invariance) . These correlations can be used to determine the universality 
classes of such models. The distribution of the absolute distance x between successive minima is a good example for 
the spatial correlation and defined by Pj U mp(x) oc , where 7r ~ 2.92 for the isotropic BS model and 7r ~ 2.57 for the 
anisotropic BS model on a square lattice. The temporal correlations Pfi rs t{t), the distribution of first return times, 
and P a u(t), the distributions of all return times scale as Pfi ra t(t) oc t~ T f irst and P a ii(t) oc t~ Tal1 , where Tf irst ~ 1.24 
and T a u — 0.71 for 2-dimensional isotropic BS model; Tf irst ~ 1.32 and r a u ~ 0.85 for 2-dimensional anisotropic BS 
model. These different values of n, Tfi rst and T a ii suggest that the isotropic and anisotropic BS models belong to 
different universality classes. 

The SOC feature of the BS model is revealed in its ability to naturally evolve towards a scale invariant stationary 
state where the correlation length is infinite and an initial local perturbation might lead to a global effect. Therefore, 
the sensitivity to the initial conditions in the BS model is a crucial matter. In order to study the effects of an initial 
perturbation one can borrow the damage spreading technique from the dynamical systems theory, which has been 
used in the literature to investigate the propagation of local perturbations in 1-dimensional @, 0, H, 0, E3, EH E3] and 
2-dimensional [5J BS models. The algorithm of this technique for a 2-dimensional BS model considered in this work 
can be described as follows: (i) once the stationary state has been achieved, consider the configuration as replica 1 
denoted by f} • , (ii) produce an identical copy of j\ ■ and introduce a small damage in this copy by interchanging the 
site with minimum fitness with a randomly chosen site (denote the new replica as ff •), (iii) let both replicas evolve in 
time using always the same set of random numbers. One can then define the Hamming distance between two replicas 




where (• ■ • } stays for the configurational averages over various realizations. Measuring the evolution of the discrepancy 
between two initially close configurations under the same external noise shows that Eq. (TT]) exhibits an initial power- 
law divergence proportional to t a , where a ~ 0.48 and a ~ 0.53 for 1-dimensional isotropic and anisotropic BS 
models [l(|, and a ~ 0.83 and a ~ 0.91 for 2-dimensional isotropic and anisotropic BS models followed by a 
finite size dependent saturation regime. These scaling exponent values of the Hamming distance suggest that 1- and 
2-dimensional isotropic and anisotropic BS models are weakly sensitive to the initial conditions (i.e., a > 0). 

The above mentioned Hamming distance measurement captures the short-time dynamics of the model very well. 
However, Tirnakli and Lyra have recently introduced a new Hamming distance measurement technique to identify 
both the short-time and long-time dynamics of the model [ill ] . Their Hamming distance definition on a chain has the 
following form, 

0i(*) = (iEi#-&A ( 2 ) 

where L is the system size and j — 0, 1, • • • , L — 1. Namely, Eq. says that the measure of the Hamming 
distance D{t) between two configurations is the smallest among the L possible values of Dj(t). It should be noticed 
that for j = (i.e., Do{t)) Eq. ([2]) turns to the conventional Hamming distance measurement. This new definition 
of the Hamming distance has been used to study the convergence dynamics of two independent configurations of 
1-dimensional BS model [l3l |. 

Our task will be to study the convergence dynamics of two uncorrelated configurations of the BS model by gener- 
alizing the recently introduced Hamming distance measurement to 2-dimensional case. The numerical procedure and 
the generalization of the new Hamming distance measure to 2-dimensional BS model are given in Section [TTJ The 
short-time and long-time dynamics of distinct configurations of 2-dimensional isotropic and anisotropic BS models 
are studied under the same external noise influence. Once the model and the numerical procedure are introduced 
then the statistical properties of the convergence dynamics of two uncorrelated configurations of the isotropic and 
anisotropic BS models constructed on a square lattice are discussed through the jump activity rate and the waiting 
time distributions in Section [IIII A summary of the results in Section ITVl concludes the paper. 
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FIG. 1: The short-time evolution of the normalized Hamming distance between an equilibrium configuration and its slightly 
modified replica for three different system sizes of the isotropic (left) and anisotropic (right) BS models. After the initial 
power-law growth D(t) oc t a where a = 0.83 ± 0.03 for the isotropic case and a = 0.91 ± 0.03 for the anisotropic case [jj, 
it decreases as both configurations converges to the same sequence of random numbers. The number of realizations used in 
calculations is 200 for N = 50 x 50 and N = 100 x 100, and 100 for N = 200 x 200. 



II. NUMERICAL PROCEDURE AND TIME EVOLUTION OF HAMMING DISTANCE 



To investigate the convergence dynamics of two uncorrelated random configurations of 2-dimensional isotropic 
and anisotropic BS models we implement the BS algorithm on a square lattice with periodic boundary conditions. 
N = L x L fitness values fij (random numbers taken from a uniform distribution in the interval [0, 1]) are initially 
assigned to all sites of a square lattice of an edge size L. As mentioned in Section HI the system evolution is standard 
re-assignment of the site with the minimum fitness and its first neighbors. In our simulations, the transient times of 
different lattice sizes for the system to achieve the statistically stationary state are attained by observing the slope 
value of the time evolution of the Hamming distance until it reaches to a fixed value as suggested in Ref . [5] . The 
chosen transient times are 4 x 10 5 , 1 x 10 6 and 2.4 x 10 6 for N = 50 x 50, N = 100 x 100 and N — 200 x 200, 
respectively. 

Once the configuration achieves the statistically stationary state we implement the standard damage spreading 
algorithm given in Section |H After a characteristic time both configurations flj and /? ■ will be composed by the 
same sequence of random numbers just shifted by a random distance. That is, // j and ff j are now indistinguishable 
and can be considered as identical replicas. 

For such an implementation of the Bak-Sneppen model on a square lattice the Hamming distance measure D(t) 
defined at each time t as the smallest among the L x L possible values of Dk,e(t) (k = 0, 1, • • • , L — 1 and £ — 
0, 1, • • • , L - 1) given by 

D k ,m = (^ E iflj-fi+kj+tij- (3) 

Notice that -Do.o(^) corresponds to the measure used in previous damage spreading studies of 2-dimensional BS model 
with the same scaling exponents [5j . 

In Fig. |T]) our results for the short-time evolution of the normalized Hamming distance D(t)/D(l) for three different 
system sizes N = 50 x 50, = 100 x 100 and N = 200 x 200 are shown in the case of 2-dimensional isotropic and 
anisotropic BS models. The initial power-law regime where D(t) oc t a extends for longer periods as the configuration 
size increases. The corresponding scaling exponents are a ~ 0.83 for the isotropic case and a ~ 0.91 for the anisotropic 
case and these a values are in good agreement with the previously reported values Q. After the initial power-law 
regime ends up at a characteristic time which depends on the system size, Hamming distance reaches a maximum and 
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FIG. 2: The long-time evolution of the Hamming distance between two uncorrelated equilibrium configurations for three different 
system sizes. After the initial transient is achieved D{t) exhibits a non-trivial decay followed by an exponential relaxation due 
to finite size effect. Additionally, from the figures one can notice that after the initial transient the Hamming distance is close 
to its initial value, D(0) — 0.224 and D(0) = 0.187 for the isotropic (left) and anisotropic (right) cases, respectively. The 
number of realizations used in calculations is 200 for N = 50 x 50 and N = 100 x 100, and 100 for N = 200 x 200. 



starts to decrease. The topological behavior of the short-time evolution of the Hamming distance in 2-dimensional 
BS model is similar to the behavior obtained in 1-dimensional BS model [Tlj |. 

In order to study the long-time dynamical regime of the Hamming distance evolution one needs to use a slightly 
different approach than the one mentioned above for the short-time dynamical regime. The new approach can be 
described as follows. Consider two initially uncorrelated random configurations of the BS model on a square lattice 
evolving under distinct noises. These configurations eventually reach uncorrelated statistically stationary states. For 
very large chains, it is well known that the Hamming distance between such configurations is 1/3 of the width of the 
fitness distribution in equilibrium This fact remains true for the configurations described on a square lattice. Then 
the initial Hamming distance between 2-dimensional two configurations of the isotropic and anisotropic BS models is 
given by 

l iff °- 224 (isotropic BS) 

m= E \flj-fU = VH (4) 

i,j=i d [0.187 (anisotropic BS). 

Once two uncorrelated configurations at their statistically stationary states are obtained we let both configurations 
evolve in time under the influence of the same external noise. That is, the minimum fitness and its nearest neighbors 
of each configuration are replaced by the same fitness values chosen from the same uniform random distribution 
in the interval [0, 1]. As time increases, both configurations converge to the same sequence of fitness values while 
displaced in space. To investigate the convergence of both configurations, the Hamming distance measure defined as 
the minimum among all displaced distances Dk,e{t) (see Eq. (|3"])) is used. As emphasized before, making use of this 
Hamming distance measure enables one to observe the long-time evolution of the minimal Hamming distance. Such 
an observation is shown in Fig. ([2]). As it is the case for 1-dimensional case, the evolution of the minimum Hamming 
distance follows the solution of a generalized nonlinear equation, firstly proposed in the context of non-extensive 
thermostatistics [Til [TBI . Il6j |. which includes in a single expression the transient, pre-asymptotic and asymptotic 
exponential relaxation regimes [HI, From Fig. @ it is seen that the pre-asymptotic regime is very close to a slow 
logarithmic decay which is associated with the uncorrelated nature of the local variables in the statistically stationary 
states. Comparing our results of 2-dimensional BS model with 1-dimensional case [ill, EH reveals that the long-time 
evolution of the minimum Hamming distance measure of both versions of the BS model exhibit topologically similar 
behavior. 



5 




FIG. 3: (color online) The time evolution of the activity rate dN/dt for three different system sizes of the isotropic (up) and 
anisotropic (down) BS models. During the initial regime, which depends on the system size the averaged activity rate is roughly 
constant crossing over to a regime of continuously decreasing activity. The initial regime scales as 1/N while the regime with 
the decreasing activity scales as 1/t (see insets). The number of realizations used in calculations is 200 for N = 50 x 50 and 
N = 75 x 75, and 150 for N = 100 x 100. 



III. ACTIVITY RATE AND WAITING TIME DISTRIBUTION 



The smallest Hamming distance among all possible displacements determined by the fitness sites {k,£} given 
in Eq. ([3]) is computed by comparing the state of two replicas introduced in the previous section. The optimal 
displacement determined by the fitness sites {k*,£*} evolve in time intermittently jumping between a few positions 
which provide almost the same Hamming distance. As it happens in 1-dimensional BS model [l3j |. after a long time 
evolution the optimal displacement between two 2-dimensional replicas stops jumping (i.e., activity) once the replicas 
becomes identical without the need of any additional displacement. 

To characterize this process, we have computed the activity rate dN(t)/dt, where N(t) is the total number of jumps 
the optimal displacement has made until time t for 2-dimensional isotropic and anisotropic BS models. The number 
of realizations used in calculations is 200 for N = 50 x 50 and N = 75 x 75, and 150 for N — 100 x 100. At this point 
it should be emphasized that the computational difficulties (i.e., elapsed time to get the results) led us to work with 
100 x 100 lattice as the largest lattice size, because the time period needed to get the most realistic results for the 
activity rate and waiting time distributions increases rapidly as the lattice size increases. 

The physical interpretation of the picture given on Fig. ([3]) is as follows: the first regime where the time evolution 
of the activity rate is roughly constant giving a plateau corresponds to the initial relaxation towards the most 
predominant displacements providing the set of smaller Hamming distances. The activity rate scales as 1/N, namely, 

dm -\f(D- <» 



dt N \N r 

During this period the two replicas stay uncorrelated and the fitness updates occur at uncorrelated random sites. 
This means that one can expect jumps to take place at a constant rate. These jumps match the minimal Hamming 
distance condition and evolve just a few displacements for which the Hamming distance is close to the smallest one. 
The second regime where the jump activity rate decreases is a consequence of the decrease in the minimal Hamming 
distance and reflects the locking of the two replicas at a predominant displacement. Configurational fluctuations of 
the time required to achieve the final displacement gives rise to the slowly decaying average activity rate characterized 
by a power-law on the form dN(t)/dt oc t . The size and time scaling of the jump activity rate is demonstrated in 
the insets of Fig. 

Comparing our results in 2-dimensional BS model with the results previously obtained in 1-dimensional BS model 
[l3| reveals that the topology of the activity rate evolution in time is almost the same for both dimensions while the 
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FIG. 4: (color online) The time evolution of the waiting times distribution for three different system sizes of the isotropic (left) 
and anisotropic (right) BS models. For both versions of the BS model it scales as f (At) oc At -2 . The number of realizations 
used in calculations is 200 for N = 50 x 50 and TV = 75 x 75, and 150 for N = 100 x 100. 



finite-size scalings are different. In our case the activity rate scales as 1/N while it scales as 1/N 2 in 1-dimensional 
case [Hj|. 

Additionally, the waiting time distribution P(Ai), i.e., the distribution of time intervals between two successive 
jumps is calculated to investigate the intcrmittence of the displacement jumps. Our results for P(At) is shown 
in Fig. ([!]). As it is the case for 1-dimensional BS model, P (At) follows a well denned power-law scaling with 
P(At) oc At -2 for both, isotropic and anisotropic cases. This means that the jumping process has no characteristic 
time scaling besides that due to finite-size effect. 



IV. CONCLUSION 



We have studied the convergence dynamics of two uncorrelated random configurations of the isotropic and 
anisotropic BS models on a square lattice by implementing the damage spreading technique. The chosen measurement 
way of the Hamming distance, which is defined as the smallest among N possible values of Dj(t) at each time step t 
enabled us to capture both, the short-time and long-time evolution dynamics of the Hamming distance. 

It has been shown that the short-time evolution of the Hamming distance exhibits an initial power-law regime 
scaling as D(t) oc t Q , where a ~ 0.83 for the isotropic and a ~ 0.91 for the anisotropic cases. These a values 
are in full agreement with those previously obtained by using the conventional Hamming distance measure. After a 
characteristic time period, which depends on the system size, the initial power-law growth of the Hamming distance 
reaches a maximum and then follows a decreasing regime. In order to study the long-time dynamics of the Hamming 
distance we have considered two initially uncorrelated random configurations of the BS model on a square lattice. 
Our extensive simulation results have shown that the initial Hamming distance D(Q) attains the value of the 1/3 of 
the width of the fitness distribution as it is the case in 1-dimensional case, i.e., D(Q) = 0.224 for the isotropic and 
D(0) = 0.187 for the anisotropic BS models characterized on a square lattice. 

Considering the algorithm based on the definition of the minimal Hamming distance between all possible relative 
displacements of two uncorrelated configurations we have reported that the optimal displacement evolves in time 
intermittently. It has been seen that in the short time period where the two replicas stay uncorrelated the activity 
rate scales as 1/N following a decreasing regime in the long time period where the activity rate has a power law, 
dN(t)/dt oc t . Moreover, the waiting time distribution was computed and shown to follow a power-law decay 
P(At) oc At~ 2 . The scaling laws obtained for the time evolution of the activity rate and the waiting time distribution 
reflect the absence of typical spatio-temporal scales in the mechanism of displacement jumps used to quantify the 
convergence dynamics. 
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The time evolution of the activity rate and the waiting time distribution for the isotropic and anisotropic versions 
of 2-dimensional BS models reveal that they both evolve in time almost the same suggesting that choosing a preferred 
direction to propagate the avalanches does not have considerable influence on the convergence dynamics of two 
uncorrelated configurations of 2-dimensional BS model. 

A difference between the BS models defined on a chain and 2-dimcnsional lattice appears on the scaling factor 
defined in Eq. ([5|). As it is seen, this factor is 1/N for 2-dimensional BS model while it is l/N 2 in 1-dimensional case. 
Although, the obtained scaling factors in 1- and 2-dimensional cases suggest a relation with the number of fitness in 
the lattice, a generalization of this relation to the higher dimensional cases needs further investigation. 
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